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CHAPTER 1 


INTRODUCTION 

1.1 Liquid Jet Instability 

The first theoretical study of the capillary (or surface tension) 
Instability of a jet of liquid was carried out by Rayleigh (1878). He 
analyzed the linear stability of an Invlscld, Incompressible, Infinitely 
long cylindrical jet (with plug flow velocity profile) to a spatially 
harmonic disturbance of the jet surface. Only axlsymmetrlc modes were 
found to be unstable and the wavelength of the most unstable (temporally 
growing) mode was determined to be 4.51 . . . times the undisturbed jet 

i - 

diameter. The cutoff wavelength for unstable disturbances, below which 
the disturbance will decay In time, was found to be equal to the undis- 
turbed jet circumference. 

Later, Rayleigh (1892) Included the effects of viscosity on the jet 
stability. He found that Increasing the viscosity had the effect of 
Increasing the wavelength of the most unstable mode. Since these stud- 
ies, numerous Investigators have studied the jet breakup process. 

Most of the analyses treated the temporal Instability problem. Lee 
(1974) used a one-dimensional, Invlscld model to study the jet breakup 
process. Proceeding on the assumption that the axial velocity depends 
only on the axial coordinate he arrived at a simplified set of equations 
for momentum and mass conservation. The linearized equations thus 
obtained were solved In closed form while numerical solutions were 
obtained for the nonlinear analysis. Among the results obtained were the 
time for the jet to breakup (In both the linear and nonlinear cases) and 
the volume ratio of the satellite drops to the total drops (in the non- 
linear analysis). In the range of parameters studied satellite drops 
were always found to occur. 
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While the temporal theory has received the most attention the actual 
jet breakup usually takes place as a result of spatially growing distur- 
bances. This was recognized by Keller, Rublnow and Tu (1973) who, after 
making a transformation to a moving jet, solved the dispersion relation 
obtained by Rayleigh with real frequency while allowing the wavenumber 
to be complex. In this way spatially growing waves are obtained. 

It was found by Keller et al. that the growth rates for the temporal 
and spatial Instabilities are the same only In the Infinite Weber number 
limit. Hence the results from the temporal theory can be applied as an 
approximation to a real jet only for large Weber numbers. 

Pimbley (1976), used the one-dimensional model of Lee to treat a 
boundary value problem where the unstable disturbances grow with distance 
from the jet exit. He found that Lee's solution (after a transformation 
appropriate to a moving jet) was the Infinite Weber number limit of the 
solution to his boundary value problem. 

Experimental work has necessarily been concerned with spatially 
growing disturbances and the jet Is seen to break Into drops at some 
distance from the nozzle rather than after some period of time. The 
results obtained for the growth rates and jet breakup lengths have gener- 
ally been converted for comparison with Rayleigh's (temporal) theory 
using the jet velocity. 

Crane, Birch, and McCormack (1964) wefe $ble to obtain uniformly 
spaced and sized drops over a range of disturbance wavelengths which 
corresponded to the range predicted by Rayleigh's linear theory. The 
growth rates of disturbances of different wavelengths were determined and 
found to depart appreciably from the theoretical predictions of Rayleigh. 
The finite amplitude of the Induced disturbances and consequent nonlinear 
effects were the reasons given for the differences. In support of this 


2 



a plot was given which shows that the growth rates In the experiment 
breaks away from the linear curve a short distance from the nozzle exit. 

Donnelly and Glaberson (1966) performed experiments on liquid jets 
of water and glycerine-water solutions In air. Disturbances were Induced 
on the jet by a nearby loudspeaker. Very clear stoboscoplc flash photo- 
graphs were taken of the jet breakup process showing the formation of the 
droplets. From these pictures the dispersion curve for the disturbances 
was determined. The spatial growth of the disturbances, as shown by the 
photographs, was converted to a temporal growth using the frequency of 
the forcing and the wavelength of the Induced disturbance for comparison 
with Chandrasekhar's (1961) linear analysis (which was based on Ray- 
leigh's analysis). The results were found to compare quite favorably. 

In light of the studies mentioned above this can probably be traced to 
the large value of the jet velocity compared with the capillary velocity 
(l.e., large Weber number) and the fact that a linear growth rate was 
found to persist to within one wavelength of the point at which drops 
begin to break off. 

Using the glycerine-water solution the effect of viscosity on the 
stability of the jet was examined. It was found that Increasing the 
viscosity had the effect of decreasing the growth rate, all other things 
being equal. These results agreed with Chandrasekhar's to within 20 
percent. Small variations In the temperature af the jet, causing the 
viscosity to vary, was cited as one reason for the discrepancy. 

A detailed review of the work done on liquid jet Instability and a 
discussion of Its applications to nozzle design has been made by McCarthy 
and Molloy (1974). Bogy (1979) has given an account of some more recent 
work concerning nonlinear effects on the jet breakup; In particular the 
formation of satellite drops. 
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1.2 Technological Applications 

A number of technologies have evolved which make use of (and helped 
to motivate) the scientific studies touched on In the previous section. 
Perhaps the best known of these Is the Ink jet printer. 

In such a printer the drops which are formed due to the jet Insta- 
bility are given a charge at the point where they break off from the jet. 
The charged drops pass through an electric field which causes them to be 
deflected from their original paths. The deflected drops then strike the 
paper and form characters. Details of the workings of such a system have 
been given by Kuhn and Meyers (1979). 

It Is easy to see that for such a system to perform properly the jet 
breakup and drop formation must be carefully controlled. A particular 
problem Is the appearance of the small satellite drops which are observed 
to occur between the main drops. For this reason a significant effort 
has been made to understand this (nonlinear) phenomenon. 

A more recent Idea for exploiting the jet breakup process Is In Its 
use as a heat transfer device, the so-called liquid droplet radiator. 

Here a heated jet of fluid Is shot out of a nozzle Into a collector a 
certain distance downstream. As the jet travels between the nozzle and 
the collector (and breaks Into drops) It radiates energy Into the space 
around It. The fluid Is then caught In the collector and circulated back 
through the system. The aim Is to develop this type of system for use 
in space on satellites or spacecraft. 

Some advantages to Its use In space over more conventional systems 
are the large surface area to volume ratio of the drops, the relative 
ease of transporting Into space a container of fluid over a network of 
pipes and Its greater chances of surviving the Impact of small particles 
floating through space. 
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A study of the feasibility of the liquid droplet radiator was made 
by Mattlck and Hertzberg (1982). While their results showed promise, 
more work Is needed before a practical system can be made. 

We have cited two examples where the liquid capillary jet breakup 
process can be put to use. These examples Illustrate the Importance of 
understanding and being able to control the jet breakup. 

1 .3 Receptivity 

It has been found In a number of problems In which a flow Is subject 
to Incident disturbances (Jones and Morgan (1972), Crlghton and Lepplng- 
ton (1974), and Goldstein (1981)) that elgensolutlons which grow exponen- 
tially downstream must be added to a particular solution In order to 
satisfy either edge requirements or causality (or both). The causality 
condition requires that the flow not respond to a disturbance before that 
disturbance Is Imposed. The edge condition which Is usually specified 
Is the Kutta condition that the velocity and pressure be finite at any 
edges or other singular points. The appropriateness of the Kutta condi- 
tion has been examined by Including viscous effects near a singular 
point. A review of work concerned with the Kutta condition has been 
given by Crlghton (1985). 

In one such study, Rlenstra (1981) using some results from triple 
deck theory obtained by Brown and Daniels (1975), showed for a semi- 
infinite vortex sheet with a plane harmonic wave Incident on It, that the 
solution satisfying the Kutta condition Is obtained as the leading order 
term In the asymptotic expansion of the outer solution for the corre- 
sponding viscous problem. Hence the Imposition of the Kutta condition 
In the Invlscld solution Is consistent with the detailed structure of 
the viscous flow near the edge. 
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The eigensolutions which are introduced into the solutions in this 
way contain terms which involve instability waves. The Kutta condition 
then provides a mechanism by which the external disturbances can couple 
to instability waves. A similar coupling occurs when causality is 
Imposed. Since specific constant multiples of the eigensolution are 
needed to satisfy each of these conditions the amplitude of the insta- 
bility waves is thereby specified. 

The triggering of instability waves by external disturbances is the 

* 

so-called receptivity problem (Morkovin (1969)). The aim of the study 
of receptivity is to determine the effectiveness of particular distur- 
bances in exciting instability waves in a flow. One measure of this 
efficiency is the amplitude of the instability wave produced by the dis- 
turbance per unit amplitude of the forcing. This quantity is known as 
the coupling coefficient. Since the Kutta condition determines the 
amplitude of the instability wave we say that the Kutta condition com- 
pletely specifies the receptivity problem. Similarly the causality con- 
dition can also be used to specify the receptivity problem. In some 
problems these conditions lead to the same result (Crighton and Lepping- 
ton (1974)) while in others they do not (Goldstein (1981)). We will be 
addressing this question for the problem studied here. 

1.4 Scope of the Present Work 

The process of the breakup of a liquid jet begins at the nozzle exit 
where spatially growing capillary Instabilities are excited by the 
external disturbance environment. The disturbances may be those which 
are coincidentally present in a facility or those due to a known forcing 
which is intentionally Imposed on the jet. Since only the axisymmetric 
mode is unstable the jet deformation is varicose ip nature (Rayleigh 
(1878)). At the point downstream of the exit where the amplitude of the 
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Instability waves become equal to the original radius of the jet a drop 
of liquid Is pinched off from the jet. Beyond this point the flow con- 
sists of a series of such drops possibly with a set of smaller drops 
(called satellite drops) occurlng between the main drops. A simple 
schematic of this process Is given In figure 1,1. 

In applications such as those discussed In figure 1.2 It Is neces- 
sary to be able to control the drop formation. In some cases for example 
It may be desired to produce drops of uniform size or to elllmate satel- 
lite drops. Control over the breakup process can be accomplished by 
Imposing an appropriate external disturbance on the flow. Since the jet 
breakup process begins with Instability waves the need to understand how 
an external disturbance effects the drop formation leads one to a con- 
sideration of the receptivity problem for a liquid jet. 

In the present work we consider the Invlscld, Incompressible, 
parallel flow of a liquid jet emerging Into a vacuum from a circular 
cylindrical nozzle subjected to a small, time harmonic disturbance. 
Specifically, this disturbance will take the form of a pulsating axial 
pressure gradient. Linearized equations for the fluctuations produced 
by this perturbation are solved subject to the kinematic and dynamic 
boundary conditions on the duct walls and free surface of the jet. 

We will require that this "steady-state" (time-harmonic) solution 
be causal. That Is, that It can develop as the long time response of the 
flow to a forcing started Instantaneously at some Initial time. Our goal 
Is to determine when spatially growing Instability waves can be excited 
by the external disturbance. 

It will be seen that there are an Infinite number of solutions to 
the problem outlined above. These solutions will differ In their 
behavior near the trailing edge of the duct. Out of these we will 
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choose the least singular one since this Is the one which can be expected 
to match to an "Inner” solution which takes viscosity Into account (Van 
Dyke (1964)). 

We begin In Chapter 2 with the general formulation of the problem. 
Physical variables are defined and the flow geometry Illustrated. We 
then obtain the linearized differential equations and boundary conditions 
for the fluctuating variables. 

In Chapter 3 we construct a causal solution to the forced problem 
using a method given by Briggs (1964). A brief description of the method 
Is given In the appendix. 

Under certain conditions the causal solution will Involve a term 
which grows exponentially downstream. This term represents a Rayleigh 
Instability wave. A major purpose for this study Is to determine the 
amplitude of the Instability wave relative to the amplitude of the forc- 
ing. This Is the so-called "coupling coefficient". 

At the end of Chapter 3 we write down a noncausal particular solu- 
tion to the same equations and boundary conditions. We will find that 
this solution Is singular at the trailing edge of the duct. 

In Chapter 4 we construct an elgensolutlon to the problem and use 
It to eliminate the singularity In the noncausal particular solution so 
that the resulting solution satisfies the Kutta condition. This solution 
again Involves Instability waves and we can compute the coupling 
coefficient. The solution obtained In th1$ way 1? compared with the 
causal solution. 

Finally, In Chapter 5 we discuss the various solutions and summarize 
the results obtained. 
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CHAPTER 2 


FORMULATION OF THE PROBLEM 


2.1 Formulation 

We consider the Invlscld, Incompressible, parallel flow of a liquid 
jet emerging from a semi-infinite circular duct of radius a Into an 
evacuated region. The flow geometry and coordinate system used Is shown 
In figure 2.1 . 

A steady base flow (which satisfies the Invlscld equations of motion 
for any choice of U(r) 


U - (U(r), 0) 


( 2 . 1 ) 


P = C = constant 

Is subjected to a small, time harmonic, axlsymmetrlc perturbation. 

Assuming the equations can be linearized the fluctuations due to 
this disturbance will likewise be harmonic In time and we write them as 


( u 1 , v 1 , p 1 ) = e 1o)t (u'(x,r), v'(x,r), p * ( x , r ) ) (2.2) 

where w Is the frequency of the disturbance. 

We write the Instantaneous variables as 

u^ = U +• u 1 

V T = v 1 (2.3) 

p T = p 1 + C 

The Instantaneous flow variables are governed by Euler's equation 
for an Invlscld fluid 


auT 

Bt~ 


T 

+ ( u • V) u 1 


= -vp 


(2.4) 


and the Incompressible continuity equation 

V • u T = 0 (2.5) 

The kinematic boundary condition on the solid duct wall Is 
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v T (x,a,t) = 0; 


—00 


( 2 . 6 ) 


< x < 0 
0 < t < 00 

The perturbed free surface of the jet can be described by an 
equation of the form 

r = C(x,t) + a (2.7) 

where £ Is the displacement of the surface from Its undisturbed 
position. 

On the free surface we have the kinematic boundary condition that 

particles on the surface move with the surface. That Is, 

\ • « 

vT(x,f. + a.t) . ot (2-8) 

where D/Dt Is the convective derivative 

St = ft + uT(c * a) k (2 - 9) 

In the absence of viscosity the dynamic condition at the free 
surface requires that surface tension forces balance pressure forces. 

This condition can be written as 

P T (x,C + a.t) = y(~ + (2.10) 

where y Is the surface tension and and are the principal 

radii of curvature of the surface at position x and time t (Landau 
and Llpshltz (1959) p. 231). 

Substituting equation (2.3) Into equations (2.4) to (2.6) and 
linearizing for small disturbances we get 


au' 
p [at 

- ♦ U(D “ 

+ ^ v . . 

IElL 

ax 

(2.11) 

p 

[ff + U < r > 

av'l 

ax J ’ " ar 


(2.12) 


V • 

u' = 0 


(2.13) 

V 1 

(x,a,t) = 0 

i 

8 

A 

X 

A 

0 

(2.14) 
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Since we have assumed that the disturbances are small, the displace- 
ment of the surface from its original position will be small. We can 
then expand the boundary conditions for x > 0 about r = a for c 
« 1. This gives 

v'(x,a,t) = |^ + U(a) 0 < x < co (2.15) 

for the kinematic condition and 


p'(x,a,t) = - y 



0 < X < co 


(2.16) 


for the dynamic condition when the disturbances are axisymmetric (Lamb 
(1945) p. 473). Note that in equation (2.16) we have set C = y/a which 
is the equilibrium condition when there is no surface displacement. 

We wish to determine the solution to the above problem for the par- 


ticular case of a disturbance corresponding to a time harmonic axial 

' - 1 0)t 1 

pressure gradient -P Q e where P Q is a constant. 

Boundary value problems with discontinuous boundary conditions such 
as we have here can be solved by the Wiener-Hopf method (Noble (1958), 


Roos (1969)). This method makes use of the analytic continuation of 
unilateral Fourier transforms into the complex plane. We make use of 


this method in the present work. 
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CHAPTER 3 


FORCED JET SOLUTION 

3.1 Introduction 

Since the present formulation Is linear we may obtain the solution 

to the forced jet problem by superposition. 

In particular, we first consider the flow In a doubly Infinite 

circular duct subject to an axial pressure gradient -P Q e where 
» 

P Q Is a constant. The solution to the doubly Infinite duct problem 
will satisfy the differential equations (2.11) to (2.15), the boundary 
condition on the solid duct wall (eq. (2.14)) and the condition that at 
upstream Infinity there Is only a time harmonic pulsation of the flow due 
to the Imposed pressure gradient. The two free surface conditions 
(eqs. (2.15) and (2.16)) however cannot be satisfied by this solution. 

For this reason we construct another solution to equations (2.11) to 
(2.13) which corrects the doubly Infinite duct solution for the presence 
of the free surface and which vanishes as x -» -<» so that the upstream 
boundary condition remains satisfied. 

In this way then the solution to the forced semi-infinite jet prob- 
lem Is obtained. 

3.2 Doubly Infinite Circular Duct 

We Impose a time harmonic axial pressure gradient aP D /ax = 

-P Q e 40 on the flow In a doubly Infinite circular duct and seek 
solutions of the form 


D, 4-X 
u (x,r,t) = 


*D 

u U (r) 


1 (\X-wt) 


(3.1) 


Since P D Is a function of x and t only equation (2.12) gives 
~D 

that v (r) ~ 0. In particular the kinematic boundary condition on the 
solid duct walls Is Identically satisfied. In order to satisfy 
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equation (2.11) we set X = 0 and u D (r) = IP ' Q /pu a constant. Hence 
the solution to this problem becomes 


P°(x,t) 

v°(x,t) 

u D (x,t) 


= -P'xe 
o 


•lot 


P o e 


-lot 


E 0 
IP' 

0 -lot 

_ e 

po 


(3.2) 


for all x where P Is the level of the pressure fluctuation at 

o 

x = 0 which will be determined as part of the solution. 

From equation (3.2) we can see that the boundary conditions of equa- 
tions (2.15) and (2.16) are not satisfied by this solution. We next 
construct a solution to equations (2.11) to (2. 13) which will correct 
equation (3.2) at the free surface of the jet. 

3.3 Semi-Infinite Jet 

In order to correct equations (3.2) for the presence of the free 
surface we will seek functions u(x,r,t), C(x,t) and p(x,r,t) f which 
satisfy equations (2.11) to (2.14), such that u^ + u and p^ + p 
satisfy the satisfy the boundary conditions for x > 0 (eqs. (2.15) 
and (2.16)). 

Using equation (3.2) In equations (2.57) and (2.16) we write 

v(x,a,t) = /ft * U(a) f(x,t) (3.3) 

and 

2 

p( x,a, t) = P^xe" 1o,t - P o e~ 1wt - y * ^~^C(x,t) (3.4) 

for x > 0. 

Then the boundary value problem which determines the desired 
functions Is 
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au 

at 


U(r) 


au 

ax 


> U ' ( r) v 




afi 

ax 


(3.5) 


av 

at 


+ U(r) 


av 

ax 




ar 


(3-6) 


<3 

• 

ic; 

ii 

o 

(3.7) 

v(x,a,t) = 0 -co < x < 0 

(3.8) 

u(x,r,t) -> 0 as x -> - » 

(3.8) 


together with equations (3.3) and (3.4). A boundary value problem such 
as that above with discontinuous boundary conditions can be solved by 
Fourier transforms using the method of Wiener and Hopf. 

Since for Incompressible flow the absolute level of the pressure 
fluctuations does not vanish at Infinity but Is felt for all x < x 
< a> we have subtracted It out from the pressure above so that the 
Fourier transform of p will exist In the usual sense. 

We will seek solutions of the form 

(u, v,p,C) = (u( x, r) , v(x,r), p(x,r), C( x) ) e" l6>t (3.9) 

where the functions u, v, p, and £ are sufficiently well behaved at 
Infinity so that their Fourier transforms In x exist. These are 
defined as 


u(k,r) 


1 _ 

2ir 


f u( x, r)e 


dx 


(3.10) 


etc . 

The Wlener-Hopf technique makes use of the analytic continuation of 
half range Fourier transforms Into the complex plane. These unilateral 
transforms are defined as 


u + (k,r) = 


1 _ 

2ir 


( 


u(x,r)e 


-Ikx 


dx 


(3.11) 
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and 


u (k,r) 



— Iky 
u( x, r)e ,KX 


dx 


(3.12) 


The Integral In equation (3.11) can be shown to converge uniformly In any 
closed and bounded subset of the upper half plane (Im k > 0) and hence 
represents an analytic function In the upper half plane (Roos (1969)). 
Likewise, equation (3.12) represents a function which Is analytic In the 
lower half plane. Along the real axis we have 


u(k, r) = u + (k,r) + u_(k,r) (3.13) 

Similar relations can be written down for the other variables. 

With the above definitions we Fourier transform equations (3.5) 
to (3.7) and obtain 

/V A 

pi ( U( r ) k - w)u(k,r) + pU’(r)v(k,r) = -1kp(k,r) (3.14) 


d p ( k ijr). 

P 1(U(r)k - <o)v(k,r) = - Hr 


(3.15) 


- v(k.r) dv( k. r) 

1ku(k,r) + r + dp =0 (3.16) 

Applying the half range transforms to the boundary conditions (eqs. 

(3.3), (3.4), and (3.8)) and noting that £(0,t) = 0 we can write 

v + (k,a) = 0 (3.17) 


v (k,a) = 1(U(a)k - w)£ 


(3.18) 


P (k,a) 



P 

o 


2ir1(k - 1c*) 


$ ‘ ^ 


(3.19) 


1c*X 

In equation (3.19) we have added a small amount of damping e 
for x > o to the forcing terms so that the Fourier transform will 
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exist. The damping factor € * will be set equal to zero at the end 
of the analysis. 

Equations (3.14) to (3.16) can be combined to arrive at Rayleigh's 

A* 

equation for v which becomes 

v(k,r) = 0 
(3.20) 

for an arbitrary base flow profile U(r). In equation (3.20) 1 

means differentiation with respect to r. 

In the following we will take 

U(r) = U^l - b 0 < b < 1 (3.21) 

where U Is a constant. Using the parameter b we can get results 
o 

for a range of profiles from plug flow to Hagen-Pol seullle flow. 

Note that for this family of profiles 


U(r) - 


O) 


r,~ 

v'(k,r) + 


v(k,r) 

r 


- k v(k,r) 


k 'AH 


1 d_ [~ U 1 ( r) ~ 
r dr r 


With equations (3.21) and (3.22), equation (3.20) becomes 


1 


1 


v " ( k , r ) + “ v 1 (k,r) - |^k + v(k,r) = 0 

Equation (3.23) Is the modified Bessel equation of order one 
Hence we can write 


(3.22) 


(3.23) 


v(k,r) = A(k) I^kr) + B(k) K^kr) (3.24) 

where 1^ and K ] are the modified Bessel functions of the first 
and second kind, respectively. 

In order for the solution to be bounded at r = 0 we must have 
B(k) = 0 so that 

v( k, r) = A(k) I ] ( kr) (3.25) 
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The function A(k) Is determined from the boundary conditions by 
the Wlener-Hopf technique. 

From equations (3.13), (3.17), and (3.25) we can write 

A A . 

v(k,a) = A(k) I-! (ka) = v_(k,a) (3.26) 

and 

A A A 

v‘(k,a) = A( k ) k I -j (ka) = v;(k,a) + v_^(k,a) (3.27) 


for k along the real axis. (In eq. (3.27) I-j(ka) = d/d(ka) 
I-j(ka)) . 

Eliminating A(k) from equations (3.26) and (3.27) gives 


v (k,a) = 


I -| ( ka ) 
kl|( ka) 


i;(k, 


a) + v 1 (k,a) 


(3.28) 


The boundary conditions (eqs. (3.18) and (3.19)) can be used to 
eliminate C to get 


. ,.2 1 

It k - 


P - (k,3) = 2ir( k -°ie*) 2 + 2 ^ k ' + t“- kU( *) ] Vjk 


.a) ♦ 


(3.29) 


Writing equation (3.14) at r = a using equation (3.29) and 
the fact that p(k,a) = p + (k,a) + p (k,a) along the real axis gives 


1p[U(a)k - a)] u(k,a) + pU'(a) v(k,a) = -Ik 


P + (k,a) - 


P' 

o 


2ir( k - 1 e* ) 


It k 


2 1 


V (k,a) + 


2ir1(k - 1 e* ) [co - kU(a) ] 2ir 

Writing equation (3.16) at r = a and using equations (3.26) and 
(3.28) we get 


(3.30) 


Ik 


kl ' ( ka) 1 

u(k,a) = -v_(k.a) [- ♦ -i^kJTJ 


(3.31) 
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Substituting equations (3.26) and (3.31) Into equation (3.30) and 
rearranging gives 

, . . 1k 2 P' k 2 (k - U*)P 

x<k)(k - U*) v_( k,a) + Ik^k - 1e*r p + (k,a) = -^T" 2 - ^ 

_ lYk 2 (k - 1e*) 2 C (0) (332) 


where 


x(k> = 


1 


[« - kU(a) ] 


p[o) - kU( a) ] 


1 kl^(ka) 
— + 


La I -j ( ka ) J 


+ pU'(a)k[o> - kU( a ) ] 


* rk 2 L - k 2> 

\a j 


(3.33) 


The key to the success of the Wlener-Hopf technique lies In find- 
ing a factorization for x(k) such that 


x(k) = 


x + (k) 

roo 


(3.34) 


where x+ Is analytic for Im k > 0 and x_ Is analytic for 
Im k < 0. 

Formally the Wlener-Hopf equation can be written as 

(k - 1c*) 2 v_(k,a) 1k 2 (k - 1e*) 2 P + (k,a) 

+ 


x (k) 


x + (k) 


1k 2 P' k 2 (k - 1.e*)P rt 1yk 2 (k - 1e*) 2 C'(0) 


2ir 


2ir 


2ir 


J x + (k) 


(3.36) 


Due to the nature of the kernal function x(k), specifically that 
It Is meromorphlc, the factorization (eq. (3.34)) can be performed In a 
simple way by making use of the Welerstrass factorization theorem (Roos 
(1969)). As mentioned In the Introduction, the solution we are seeking 
Is the least singular causal solution. In view of this we will construct 
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our factorizations so as to be consistent with the requirements of 
causality. 

3.4 Factorization of the Kernal Function and Inhomogeneous Terms 

3.4.1 Causality 

A general method for obtaining a causal response of a system to an 
Impulsively started forcing has been given by Briggs (1964). Our concern 
here Is that our "steady state" (time harmonic) solution be the long-time 
response of the flow which develops In a causal way from a forcing which 
was Initiated at some time, say t = 0. The causality condition Is that 
there Is no response before the forcing Is "turned on", that Is for t < 0. 

So as not to disrupt the discussion of the analysis a description 
of the main points of Briggs' method Is deferred to an appendix. At this 
time we merely point out that a causal solution can be obtained by solv- 
ing the problem with the frequency « having a large positive Imagi- 
nary part. The (causal) solution for real o> Is then obtained by 
analytic continuation. 

The requirement that Im be large has no effect on the algebra 
leading to equation (3.32). The factorization of the kernal function, 
however. Is effected In a very Important way. We now proceed to deter- 
mine this factorization. 

3.4.2 The Welerstrass Factorization Formula 

The Welerstrass factorization formula allows us to write an entire 
function with simple zeros as an Infinite product (Roos (1969) p. 174) 


E(Z) 


Z 

= E(o) e 


d/dZ{log[E( Z) ] } 


Z=0 

n= 



(3.37) 


where the a n are the (simple and nonzero) zeros of E(Z) provided 
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that 


00 

Eft 

n=l \ , 


< 00 


(3.38) 


for any R > 0 (Conway (1978) p. 170). 

If we represent the zeros of E In the upper half plane by 

a U and those In the lower half plane by aj" we can rewrite 
n n 


equation (3.37) as 


E(Z) = E(o) e 


Z[log[E(Z)]] Z=Q 


n=l 


•7/ U 

7 \ Z/a r, 

1 - — ^e n 
a n 


n= 


We can then define functions E + and E_ such that 


7 \ z/a n 

1 - -f ie n 

a n 

(3.39) 


E(Z) = 


E t (Z) 
E. (Z) 


(3.40) 


where 


E + (Z) = e 


4>( Z) 


n = 


Z/a 


1 - 


(3.41) 


and 


E (Z) = 


♦ (z) . - z P°g[E(z )]] z=0 


Z/a 


E(o) 


(3.42) 


n=l 


In equations (3.41) and (3.42) E + (Z) Is analytic and nonzero In the 
upper half plane and E_(Z) Is analytic and nonzero In the lower half 
plane. d(Z) Is an entire function chosen so that E + Is algebraic 
at Infinity (Noble (1985) p. 15). 
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To arrive at the desired factorization for our Kernel function we 
first write x(k) In dimensionless form and group terms as 

1 

x 2 n-x 2 ) V** ~ ~ * u o 

+ — * ' - — 2bx[« - x(l - b) ] — - (3.43) 

x x a 

or 




where 


x(x) 


1 F l<*> ^ 

[Q - x(l - b) ] F 2 (x) a 


(3.44) 


~ ~ 2 - x 2 ) *i^ x ) 

F,(x) = + I ( x ) [Q - x(l - b ) y + -~- 

\ x / B x 

Ii(x) 

- 2bx[Q - x(l - b) ] - L 7 ~ (3.45) 

x 


and 


I. (x) 

F ? (x) = -V"' (3-46) 

x 


In equations (3.43) to (3.46) we have defined dimensionless 

2 2 

quantities as Q = ua/U , x = ka and 6 = pall /y where Q Is 

o o 

2 

the Strouhal number and 3 Is the Weber number and we have used equa- 


tion (3.21) for U(r). Since the Imaginary part of w (and hence of Q) 


Is positive the term ^/[Q. - x(l - b)] Is analytic In the lower half 


plane. The entire functions and can be factored using the 

Welerstrass formula. 
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Consider first the function F^. Let 8 n be the zeros of 
(which are the zeros of 1^) In the upper half plane. Since F^ Is an 
even function of x the zeros In the lower half plane are simply -S n * 
All of the zeros of F^ are pure Imaginary. If x = 1 <j Is pure 
Imaginary we can rewrite ?2 as (Abramowltz and Stegun (1964)). 


MW) J,<«) 
f?' 1 -) > -T«- ' « 


(3.47) 


As cr -* ® . we have 


y*) 


COS a - 


3ir 


(3.48) 


Since the zeros of the cosine grow like n as n -> «>> the 
relation (eq. (3.38)) Is satisfied and equation (3.37) can be applied 
We can easily find that 


F 2<°> * 2 


(3.49) 


and 


d_ 

dx 


log[F 2 (x)] 


= 0 


x=0 


By equation (3.39) then we can write 


(3.50) 


F 2 (x) - 2 


n=l 


1 - 



1 * 


x 


(3.51) 


The function F^(x) can be handled In a similar way. Let 
be the zeros of F^ In the upper half plane and C m those In the 
lower half plane. Asymptotically the zeros of F^ becomes just the 
zeros of 1^ . In particular they grow like m as m -> <». We can 
evaluate the constants needed In equation (3.37) as 

F-j ( 0) = Q 2 (3.52) 


22 



and 


d_ 

dx 


|log[F 1 (x)']1 


b - 2 


(3.53) 


x=0 


Then using equation (3.37) 




m= 


v w x , x[l/f tl/5 ] 

1 - p- ) f 1 - jj- | e (3.54) 

Sn / V *m 


Combining equations (3.44), (3.51), and (3.54) we can rewrite 
the kernel function as 


2n 2 e x(b-2)/n 


x(x) = 


m= 


1 X If, X ! X < 1/ V 1/ V /' U 0 


[Q - x(l - b) ] 


n= 


i - Mfi + ^-U x(1/V1/Bn) 


B 


n 


n 


(3.55) 


Now let r = -<r and define 
m m 


x t (x) - 


m=l 


, x \ - x/f m 
14- — \e 


“m 


<t>(x) 


n=l 


1 + — le" X/Bn 


3 


n 


(3.56) 


and 


x (x) = 


ffi - x(l - b) 1 n- 


.11 - ^le X/Bn 


6 


n 


x4>( x) 


faU \ 2Q 2 e x(b 2)/Q 



m= 


1 - ^ le^ 


(3.57) 


so that 
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x(x) = 


(3.58) 


x + (x) 
xjx) 

Equations (3.56) to (3.58) then give the desired factorization. 

In order to determine the function 4»( x ) In equations (3.56) and 
(3.57) we need to determine the asymptotic behavior of the Infinite 
products In these equations. An asymptotic expansion for the Infinite 
products, with x In the upper half plane, has been given by Noble 
(p. 128). We can apply this expansion directly once we know the 
asymptotic form of the roots 5 n and C m - 

From equation (3.48) we can find that 

S n ~ 1ir( n + n Q ) + ~ (3.59) 

as n ■+ oo where n 0 Is some Integer. 

Recall that the zeros of F-j approach those of F^ as x ■* «. 

The asymptotic form of C then Is formally the same as equation 

m 

(3.59) so we can write 

? ~ 1ir(m + m ) + t 3 - (3.60) 

s m o' 4 

as m -> °° where m Is some Integer. The difference between the 

o 

Integers m and n will be equal to the number of roots In the 
oo 

lower half plane which are left over after we Identify the set of roots 
of F-j which can be put Into a one-one correspondence with the roots 
of F ? . The roots of these two functions were determined numerically 
and It was found that each of the roots of F ? could be associated 

with a root of F^ with four roots of F^ left over. We will denote 

the number of these four roots which lie In th$ lower half plane by 1. 

The numerical value of 8. may be a function of the particular values of 
the physical parameters being considered. Now since there are 8. roots of 
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in the lower half plane beyond those which can be associated with 

roots of F 2 we must have m Q = n Q - 8. in equation (3.60). Then 

the correct asymptotic form for £ is 

n 


f n ~ W(n + n Q ) + ~ (1 - 44) (3.61) 

With equations (3.59) and (3.61) we can write the asymptotic 
expansion for x f as (Noble, p. 128) 

*+<*> ~ exp < ♦<‘x) - 1 (1 - 41) ♦ • log tj(l- 41) + l) 



If we choose 



we will have that 


(3.63) 


X + (x) ~ x l (3.64) 

as x ^ oo. 

A numerical study was made of the "left over" zeros of F^(x) over 
a range of the physical parameters In order to determine the factor a,. 
The results of that study are now described. 

3.4.3 Numerical Results for the Roots 

Of the four roots we are considering here two are real and two are 
complex conjugates when Im Q = 0. Causality considerations require 
that we determine the position of these roots when Im Q becomes 
large. The movement of the roots of F^(x) as Im Q Is Increased from 
zero was studied numerically for a number of combinations of Strouhal 
number, Weber number and mean velocity profile. 
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For each mean profile there Is a range of Weber number and Strouhal 
number for which the four roots move as shown In figure 3.1. This sketch 
shows the movement of the roots as Im Q Is varied. The Important 
thing to note about this figure Is the movement of the complex root which 
lies In the lower half plane when Im ft = 0. For fixed profile and 
Weber number there Is a range of Re (1.e. p the Strouhal number) for 
which this root crosses the real axis and moves Into the upper half plane 
as Im Q Is Increased as shown In figure 3.1> This crossing of the 
real axis of a complex root Is the criterion for the existence of spa- 
tially growing waves In the causal solution (Briggs (1964)). The range 
of Re Q for which this crossing occurs then corresponds to the 
Strouhal number range over which the flow supports Instability waves. 
Outside of this range the flow Is stable. For a given mean profile and 
Weber number the range of Strouhal number over which the flow supports 
Instability waves can be determined by finding the range of Re (Q) 
over which this root crosses the real axis. The locus of this root as a 
function of (complex) Strouhal number was determined for a number of 
combinations of mean profile and Weber number. The results of the com- 
putations are shown In figures 3.2 to 3.19. 

For the plug flow profile (b = 0) the cutoff Strouhal number (above 
which the flow is stable) was found to be equal to one Independent of the 
value of the Weber number In the range computed. Results obtained for 
other profiles (b * 0) Indicate that for fixed b the range of Strouhal 
number over which Instability waves exist decreases with Increasing Weber 
number (see figs. 3.9 to 3.14 and 3.15 to 3.19). Additionally, for a 
given Weber numer, this Strouhal number range was found to decrease 
quite rapidly with Increasing b. In fact when the profile parameter 
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was Increased beyond 0.5 no roots were found which correspond to Insta- 
bility waves for the range of Weber number reported here. 

In the course of studying the movement of these roots another sltu- 

2 

atlon was discovered. For the plug profile when 0 < 5 the roots move 

In the way illustrated In figure 3.20. Notice that In this situation 
two roots which appear In opposite half planes when Im ft » 1 swap 
their positions for Im ft = 0 when the Strouhal number Is Increased 
beyond a certain value. Using the criterion of Briggs (1964) this Indi- 
cates that the flow Is absolutely unstable for these Weber numbers. That 
Is, disturbances grow In time at every point In space. When the flow Is 
absolutely unstable It Is no longer appropriate to consider spatially 
growing instability waves. 

For the purposes of this study, since we are Interested In the gen- 
eration of spatially growing Rayleigh Instabilities, we will confine our 
attention to the parameter ranges for which the flow Is convectlvely 
unstable. In this range the roots move as In figure 3.1 so that 8. = 2. 

With the position of the roots of F^ thus determined the factor- 
ization of the kernal function Is complete. It remains only to factor 
the Inhomogeneous terms. This Is particularly simple In this case. 

From equation (3.36) we can see that the Inhomogeneous terms (those 
on the right) are by Inspection "plus" functions. We can now proceed to 
solve the Wlener-Hopf equation (eq. (3.36)). 

3.5 Solution of the Wlener-Hopf Equation 

"Finding the solution" In problems of this type consists of finding 

expressions for both v_ and p + . This Is accomplished by using 
analytic continuation arguments along with llouvllle's theorem on entire 
functions . 
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We begin by rewriting equation (3.36) as 
2 


v (k,a)(k - 1e*) 
x (k) 


1k 2 (k - U*) 2 P + (k,a) IkV 


x t (k) 


2irx + (k) 


k (k - 1c*)P 
2irx + (k) 


Wk 2 (k - 1e*)V(0) 

2irx + (k) 


(3.65) 


and recalling that this equation holds along the real k axis. 

The term on the left Is analytic In the lower half plane while the 
one on the right Is analytic In the upper half plane. Since the two are 
equal along the real axis we can consider the former to be the analytic 
continuation of the latter Into the lower half plane. In this way we 
can construct an entire function E(k) as 


v (k,a)(k - 1e*) 2 - 1k 2 (k - 1e*) 2 P (k,a) 

E(k) = — ^ = — + 


x 00 


x + (k) 


2 

k iP i 

2irx + (k) 


k (k - ic * )P o 1vk 2 ( k - 1c*) 2 f ' ( 0) 
2irx L (k) 2irx j _( k) 


(3.66) 


where E(k) Is equal to the "minus" functions for all k in the lower 
half plane and the "plus" functions for all k In the upper half plane. 

Different choices for E(k) will lead to different solutions to the 
problem. The least singular solution will correspond to the most rapid 
decay of the Fourier amplitudes (and hence of E(k)) as k -> «>. In 

order for the Fourier transform to exist the Fourier amplitudes must be 

algebraic at Infinity so that E(k) Is at most algebraic at Infinity. 
Using this fact and an extended version of the usual Llouvllle theorem 

(Noble (p. 6)) we can determine that E(k) must be a polynomial and we 

write 




E( k) = a + bk + . . . + qk 


(3.67) 
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where q Isa finite Integer. Combining equations (3.66) and (3.67) we 
have 


V_ (k,a) = 


x_(k) 

(k - U*)' 


a + bk + ck 2 + 3k 3 + 


. . + qk q j 


(3.68) 


for Im k < 0 and 


P + (k,a) 


P‘ 

o 

2ir(k - 1c*) 


IP 


2ir( k -'1c*) 


T.C( P ) 

2ir 


1x + (k) 


k 2 (k-1e*) 2 


• a + bk + ck 2 + dk 3 + 


+ qk 


] 


(3.69) 


for Im k > 0. It remains to find the coefficients a, b, etc. 

Our choice for the coefficients of the polynomial E(k) will be 
dictated by the following requirements. These are: (1) that the solu- 

A 

tlon satisfy the boundary condition at upstream Infinity (2) that P + 

Is analytic In the upper half plane, and (3) that we obtain the least 
singular solution. 

The upstream boundary condition Is that the correction to the doubly 
Infinite duct solution (which Is what we are constructing here) vanish 
as x -» -co so that we are left with only the Imposed forcing (which Is 
given by eq. (3.2)) In the complete solution. The behavior of the pres- 

/v 

sure as x -* -«> is determined by the behavior of P + as k -7 0 (Roos, 

(1969), p. 151). In order for P to vanish as x -«> the residue of 

P + at k = 0 must be zero. This Is accomplished by setting a* = b = 0. 

* 

Since we have defined the function P + to be analytic In the 
upper half plane we must retain enogh coefficients, and choose their 

values, so that the pole at k = 1c* In equation (3.69) (which lies 

In the upper half plane) will cancel out. If we expand equation (3.69) 
near k = 1e* the singular part becomes 
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1 


(k - 1c*) 


1X + d€*) 

(U*) 2 


C(1e*)^ + d ( 1 € * ) + 


. + q(1**) q l + 2 ^ 


1 


(k - 1c*) 


^ X+ (l€*) 2(1c*)x + (1«*)’ 


( 1 € * ) 


(1c*) 


• (c(1c*) 2 + d ( 1 e * ) 3 + . . . + q(1c*) q 


IP 1x.(1«*) „ n I 

+ ^ 2c1c* + , . . qq( 1 c*) q_ 

2ir (1c*) 2 V 

and the double pole at k = 1c* will vanish If we choose 


Ip. 


c = 


,3*0 


2lTX + (l€*) 

and the overall level of the pressure fluctuations as 

-1P n x. (1c*) 

P o = 


(3.70) 


(3.71) 


X + ( 1 € * ) 

the remaining contributions vanishing as c* -» 0. 

The least singular solution (near x = 0) will correspond to the 
one whose Fourier amplitude decays most rapidly as k -* «>. Hence 
the least singular solution will be obtained by choosing the order of 

the polynomial E(k) to be as low as possible while accommodating the 

* 

consistency condition (that P Is analytic In the upper half plane) 
and the boundary conditions. Since these conditions have already been 
satlslfled we will set the remaining coefficients to zero. Then the 
least singular solution becomes 


IP 


VJk.a) = 


k 2 x_(k) 


2irx + ( 1 €*) (k - 1c*) 


(3.72) 
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and 


P.(k,a) = 

2ir( k - U*) d 


1 - 


X + (k) 


X + d€*) 


P^O**) 


2ir 


(3.73) 


2ir X+ ( 1 €*) ( k - U*) 

With the boundary values of the solutions thus obtained we can now 
write down the equations for the flow variables . From equations (3.25), 
(3.26), and (3.72) we have for the upwash velocity 


v( x , r) 


ip ; 

2irx + (o) 


•'—00 


I l (kr) Ikx 
I-j (ka) e 


dk 


(3.74) 


The streamwlse velocity can then be determined from the continuity equa- 
tion (eq. (3.16)) as 


u( x, r) = 


- p ; 


2irx + ( o) 


x_( k) 
I-j ( ka) 


I-,(kr) 

kr 


+ I-j (kr) 


Ikx 

e 1KX dk 


(3.75) 


and the pressure from the momentum equation (eq. (3.14)) 
P(x,r) = 


Pi f x (k) 


kl-j ( ka) 


2irx + (o) 


JL 


/ 1-, ( kr) 

P [o) - kU( r ) ] l + I-j< kr) 


+ pU 1 ( r) 1^ ( kr) 


Ikx .. 
e dk 


(3.76) 


The solutions of equations (3.74) to (3.76) are valid only when Im 
<j> Is sufficiently large. The "steady state" solutions for o> real 
are obtained by analytically continuing these solutions to the real co 
axis. From the discussion of the last section we know that If we are In 
the range of Strouhal number for which Instability waves exist and the 
flow Is convectlvely unstable one of the poles of x ♦ we will call 
It a, moves from the upper half plane Into the lower half plane as 
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Im a -* 0. In the analytic continuation of our solutions then we must 
deform the Integration contour In the k plane (which Is originally the 
real axis) around this pole. We can evaluate the contribution of this 
pole using the residue theorem. For the upwash velocity this becomes 


v(x,r) = 


ip ; f v kr) 

2kx + (o) x - (k) I,(ka) 

•Lm 


e 1kx dk 


I I 

frr T-hj) 1im [< k - <k)] 

(0) I^aa) k ^ a 


laX 


(3.77) 


Or, making use of the factorization of the kernal function In the second 
term 


v( x, r) = 


IP 1 

o 


r 00 


2irx + (o) 


,,, V kr) Ikx 

x - (k) I^ka) e dk 


P' 

0 


X + (a) 


I^ar) e 


1aX 


(3.78) 


X + (o) I-|(aa)x‘(a) 

Equation (3.78) Is valid for o> real. 

Since a Is In the lower half plane when u Is real the second 
term on the right In equation (3.78) grows exponentially In x. This 
term represents a Rayleigh Instability wave with amplitude - p 0 x + ( a )/ 
I^aaJx 1 (a)x + (0) . We can think of this Instability wave as being triggered 
at the trailing edge of the duct by the external forcing. The ratio of 
the amplitude of the Instability wave to that of the forcing Is the so- 
called coupling coefficient which In this case becomes 


C = 


-x, (a) 


(3.79) 


0 x + (o) I-j (aa)x' (ot) 

The magnitude of the coupling coefficient Is a measure of the effective- 
ness of the external disturbance In generating Instability waves. 
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Using the factorization obtained earlier we can compute numerical 
values for the coupling coefficient over a range of the physical para- 
meters. We do this In the next section. 

Before proceeding to the numerical calculation of the coupling coef- 
ficient we wish to point out one further feature of this problem. 

We could have constructed a particular solution to equations (3.5) 
to (3.7) and boundary conditions (eqs. (3.3), (3.4) and (3.8)) without 
regard to causality. This solution would be formally the same as the one 
constructed here up to equation (3.76) . The definitions of the split 
functions x and x , however, would be different. In this case the 

T — 

correct half plane for the real roots Is determined by allowing <o to have 
a small positive Imaginary part. This small Imaginary part will leave the 
pole k = a In the lower half plane. Hence It will not appear as a pole 
of x (k). We can obtain this particular solution from the causal solu- 
tion (with Im o) large) by replacing x + (k) by (k - a)x + (k). Making this 
replacement In equation (3.72) (with e* = o) gives 

-p - 1P o 

- (k ’ a) = 2irx ( 0 ) a < k - o) *- (k) (3.80) 

for the upwash velocity. 

3 

From equation (3.33) and (3.64) we can determine that x ~ k and 
2 

x + ~ k as k -7 <» for b# 1. Then by equation (3.34) we must have that 

1 D 

x_ ~ k (for b v= 1 ) . The behavior of v p near x = 0 Is determined 
by the behavior of V_ as k -7 «> (Roos (1969) p. 151). Inserting 
the asymptotic behavior of x_ Into equation (3.80) we find that 

V p (k,a) ~ k° as k 7 
from which It follows (Llghtlll) that 

V p (x,a) ~ <5( x) as x -> 0 + 
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where S(x) Is the Dirac delta function. Hence V p has a delta- 
function like singularity at the trailing edge of the duct. In the next 
chapter we will construct an elgensolutlon to equations (2.11) to (2.16) 

P 

which has the same singularity at the edge as v . It will be shown 
that the difference between the causal solution and the singular partic- 
ular solution Is a constant multiple of this singular elgensolutlon and 
that the constant Is just equal to that required to cancel the edge 
singularity In the particular solution and hence satisfy the Kutta con- 
dition. These aspects will be discussed In more detail In the next 
chapter. 

3.6 The Coupling Coefficient 

For the numerical computations we will return to dimensionless 
variables. We make the coupling coefficient dimensionless as 


pU, 


C = 


— C = 
a o 


-x + ( °0 


pU 


(3.81) 


X+ (0)I 1 (a)x , (a) 

where x Is given by Its Infinite product representation (eq. (3.56)), 
+ 

x Is given by equation (3.43) and 1 now means d/dx. 

We can rewrite the Infinite products In equation (3.56) In a form 
more convenient for computations by using some results given by Noble 
(p. 128). Following Noble we compare the asymptotic behavior of the 
Infinite products In equation (3.56) with that of 


0(x) = j 
n= 


+ 


— * — ^ 
a^n + b.j I 


-x/a^ n 
e 



(3.82) 


where a^ and b^ are the coefficients In the asymptotic form of 
the roots, and c Is Euler's constant (c = 0.5772 . . .). 
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If we divide the Infinite products In equation (3.56) by the Infin- 
ite product representation of J and Insert the values for a^ and 
we can write 


x Z (1/lirn-l/? ) 

' n=l n 

N ~ “ e 

J (x) 


00 

- 

lirn + Y (1 - 49.) 

n = 


V 


n= 


1 + 


? n - 1*n - T 0 - «) 


kn + — ^ (1 - 48.) + x 


(3.83) 


and 


X°(x) 

J D (X) 


= e 


x L (l/lirn-l/6 ) 
n=l n 




i „ , lir 
lirn + — 
4 



5_ 

n 

n= 



n = 


1 + 


B - lirn - ^ 
n 4 


. Iir ~ 

lirn + + x 

4 


(3.84) 


where the superscripts N and D are used to distinguish between the 
Infinite products In the numerator, and denominator, respectively. 

From the gamma function representation of J we have 


,N,~ V -cx/W 
J (x) = e 


r (Lil 


"ft * ! * 


(3.85) 


J D (5) = e- cx/1ir 


r fe * ! 


From the functional equation for the gamma function 


(3.86) 


r S') . 5 

1 \4/ \4 1 


4 - 0 r (t - 1 


(3.87) 
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and 



Combining equations (3.56) (3.63) and (3.82) to (3.88) we can write a 
new expression of x + as 


00 



lir 


1 + 


C n - lirn - “ (1 - 49.) 


lir 


n=l 


lirn + (1 - 4fi.) + x 


1 + 


8 - Wn 

n 


lir 

4 


n= 


1 ir ~ 
lirn + 7 - + X 
4 


(3.89) 


The roots c n and B n approach their asymptotic values fairly 
quickly. Expressing the terms In the Infinite products so that they 
occur as differences between the roots and their asymptotic values 
reduces the number of arithmetic operations needed before the products 
can be considered converged. Not only Is the computation time reduced 
by this procedure but the error due to roundoff which would accumulate 
after many operations Is minimized. 


We can use equations (3.44) to (3.46) to write x'(£) as 


x'(a) 



F'(a) 


!-,(<*) 



(3.90) 


36 



by noting that F^a) =0. 

From equation (3.56) we see that x (o) = 1. 

■f 

Combining equations (3.81), (3.89), and (3.90) along with the fact 
that % = 2 to write the final computational formula for the coupling 
coefficient as 



We now give the results of the numerical evaluation of the coupling 
coefficient using equation (3.91). 

3.7 Numerical Results 

Numerical values for the magnitude of the coupling coefficient were 
computed for the mean velocity profiles considered In section 3.4.3 and 
a number of Weber numbers over the range of Strouhal numbers for which 
spatially growing waves were found to exist (see section 3.4.3). 

Figures 3.21 to 3.23 show the results of the computations. Each 
curve In these plots show the variation of the magnitude of the coupling 
coefficient with Strouhal number for a fixed Weber number. 
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For the plug profile (fig. 3.21) | C | Is nearly Independent of 

over the range shown while for the other mean profiles (figs. 3.21 

and 3.23) |C| Increases with Increasing Weber number. 

According to these results, for a given profile and Weber number, 

the greatest coupling generally occurs for smaller Strouhal numbers 

2 

although there are curves (b = 0.4, 3 =10 most notably) for which 

| C | reaches a minimum and then begins to Increase as the cutoff 
Strouhal number Is approached. By superposing these figures we can 
determine the effect of the mean velocity profile on |C|. When this 
Is done It can be seen that, for fixed Strouhal and Weber numbers, the 
magnitude of the coupling coefficient Increases with Increasing b. This 
result Indicates that a greater coupling could be achieved In practice 
If a longer nozzle were used since then the mean profile would be more 
fully developed. Along with this Increase of the coupling coefficient 
however comes a decrease In the growth rate of the Instability wave (see 
section 3.4.3) . 
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CHAPTER 4 


THE EIGENSOLUTION 

4.1 Introduction 

In the previous chapter, In addition to obtaining a causal solution, 
we obtained a noncausal, singular particular solution to the boundary 
value problem In section 3.6. We will show that the singularity In this 
particular solution can be removed by adding the correct multiple of an 
elgensolutlon which has the same order of singularity. The resulting non- 
singular solution Is said to satisfy the Kutta condition. It remains to 
be determined however whether or not the solution obtained using the Kutta 
condition Is causal . 

We will construct the necessary elgensolutlon In this chapter. Once 
this has been done the elgensolutlon will be combined with the particular 
solution (eq. (3.80)) In such a way as to eliminate the singularity. The 
nonsingular solution thus obtained will then be compared with the causal 
solution obtained In Chapter 3. As In the previous chapter we construct 
the elgensolutlon by superposition. 

First we consider the flow of a doubly Infinite jet of fluid In the 
absence of any duct. That Is, we find a solution to equations (2.11) to 
(2.13) subject to free surface boundary conditions for all x, -<» < x < 
a.. Since this solution (which Is a Rayleigh Instability wave) will not 
satisfy the boundary condition on the solid duct wall, a solution Is con- 
structed which cancels the normal velocity on the boundary for x < 0. 

The sum of these two solutions then Is the desired elgensolutlon. 

4.2 Doubly Infinite Jet 

In this section we seek solutions to equations (2.11) to (2.13) sub- 
ject to the free surface boundary conditions (eqs. (2.15) and (2.16)) for 
all x, -oo < x < «». We look for solutions of the form 
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where ac = a> with a and c complex. Substituting equation (4.1) 

Into the equations of motion we obtain the following equations for the 
functions 

P 1(U(r) - c) u X ( r ) + pv ! (r) = UpV) (4.2) 

P 1a(U(r) - c) 0 X ( r) = -p 1 ' (r) (4.3) 

and 


J i«u I (r) + 1 ~ £rO I ( r) ] .0 (4.4) 

These equations are formally Identical to equations (3.14) to 
(3.16) of Chapter 3 (In this case the functions are the normal mode 
solutions whereas before they were Fourier transforms). They can be 
manipulated In the same way to obtain the Rayleigh equation for v 1 as 

(U(r) - Ov^fr) + i 1 ' (r) 


- 'b [“*] *;<»«- 4 1 * y 

Substituting the family of velocity profiles we are considering 

- 1 

we can again arrive at the modified Bessel equation for v 
; IH (r) + ^ v 1 ' ( r ) - v*(r) = 0 

After applying a boundedness condition on the solution at r = 0 
obtain 



(4.5) 


(4.6) 


(4.7) 


we 


vV) = Al^ar) 


(4.8) 


where A Is an arbitrary constant. 



The boundary conditions (eqs. (2.15) and (2.16)) determine the 
values of a for which nontrivial solutions exist, that Is the 
eigenvalues. Using equation (4.1) In equations (2.15) and (2.16) we 
obtain 


v 1 < a ) = la C I (a)(U(a) - c) 


and 


(4.9) 


P*(a) Y C X (a) (a? - ^ 


(4.10) 


- T 

Eliminating C (a) from these two equations we can write 


-I 

v L a) 


1a(U(a) - c) 


*1 

-filial 


/ 2 1 
Y a - - 


(4.11) 


-I -I 

Another relation between v (a) and p (a) can be obtained by el 1m- 
- 1 

Inatlng u between equations (4.2) and (4.4) and evaluating at r = a. 

Upon substituting for the mean flow profile (eq. (4.6)) this relation 
becomes 


P^a) - =£ 1 


, / 2U b U (1 - b) - c 

'‘I, X I 0 0 

V (a) — — + 


+ (U 0 {1 - b) - c) v 1 ' (a> 


We can now combine equations (4.11) and (4.12) along with the 
solution (eq. (4.8)) to obtain the eigenvalue relation 


(4.12) 


I -j ( <*a ) 


3 

_ £a_ 


1 


[“- V 1 - b) “ ] " Y aa(aV - 1) 


u - U (b + l)a 
o v _ 

ad 


I^aa) + [o> . U o (l - b)a] I-j(aa) 


for the case b = 0 this reduces to 

1 


I-| (<*a) 


[<*> " U Q a] 


3 

_ 


Y aa(a 2 a 2 - 1) L 


G) - U a 
0 


aa 


I-j(aa) + (a) - U Q a) I-j(aa) 


(4.13) 


(4.14) 


which agrees with the result obtained by Keller et al. (1973). 
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Equations (4.8) and (4.13) give the eigenfunctions and eigenvalues 
of the normal mode solution for the doubly Infinite jet. We will now 
construct a solution which corrects the normal mode solution for the 
presence of the semi-infinite duct. 

4.3 Semi-Infinite Jet 

The solution constructed In section 4.2 has a nonzero normal veloc- 
ity on the boundary for x < 0. To correct for the existence of the 
rigid duct walls we must superpose on this solution one which cancels the 
normal velocity on the boundary for x < 0. 

Specifically, we seek functions u, v, p and c which satisfy 
equations (2.11) to (2.13), boundary conditions (eqs. (2.15) and (2.16)) 
for x > 0 and the following boundary conditions for x < 0 

v(x,a,t) = -Al^aa) e 1(oX '“ t) (4.18) 

The boundary condition eq. (4.18) ensures that the normal velocity 
at r = a Is equal and opposite to that In the solution of section 4.2. 
Applying the unilateral Fourier transform (eq. (3.11)) to the boundary 
condition (eq. (4.18)) we get 


1AI (aa) 

V k,a) " 2*(o - kj 


(4.19) 


As In Chapter 3 we will assume time harmonic solutions which can be 
Fourier transformed In x and use the Wlenqr-Hopf technique. 

The same manipulations leading to the Wlener-Hopf equation In Chap- 
ter 3 can be performed here as well. The only difference being In the 
Inhomogeneous terms In the Welner-Hopf equation (due to the different 
boundary conditions for the two problems). Without repeating the algebra 
here we simply write down the Wlener-Hopf equation for this problem as 


42 



x(k) V_(k.a) ♦ 1k 2 P t . - x(k) - r [I ( ^ k u(I)] 


lAI^aa) 
2ir(a - k) 


(4.20) 


In constructing the (singular) elgensolutlon, as In the particular 
solution, the frequency to Is allowed to have a small positive Imaginary 
part. The factorization of x(k) when Im to Is small Is related to 


the factorization when Im to Is large by 


x + (k) = x + (k) (k - a) 


(4.21) 


where x are the split functions obtained In Chapter 3 with Im to » 1. 

4 * 

Using the functions x + here we can rewrite equation (4.19) as 


VJk.a) Ik P + (k,a) 
+ 


xjk) 


x + (k) 


L(k) x 


2 2 2 1 
Ykjk' - 1,/a j . 

+ ( k) [to - kU(a) ]_ 


lAI^aa) 
2ir (a - k) 


2irx + ( k) 


(4.22) 


The last term on the right hand side of equation (4.22) Is already 
a "plus" function. The remaining Inhomogeneous terms must be factored 
by subtracting out poles. For this purpose we define 


G(k) = 


1 A I -j (aa) 


2 2 2 

L _ + YkJkj^lAO 

X_(k)(a - k) ; (k)[» - ku(a)](a - k) 


(4.23) 


and construct functions G + and G so that 

G(k) = (G f (k) + Uk)) 


(4.24) 


where G + Is analytic In the upper half k plane and G Is analytic 
In the lower half plane. By Inspecting equation (4.23) we can see that 


the desired functions are 
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G + (k) = 


~ 1AI l (aa) 1AI-, (eta) vk 2 (k 2 _ 1/a 2) 


2ir( a - k ) x_( ci ) (a - k)[u - kU(a) ]x + (k)2ir 


G_(k) = 


-1 AI^ (aa) 



(4.25) 


{yjall ((u“av) l/a^lAI^aa) 


2ir( a - k)x_(k) 2ir(a - k)x_(a) 

nte uuta) 

(- - ufiy) - ku(a)]2 " 

Substituting equations (4.25) and (4.26) Into equation (4.24) and 
rearranging we get 


(4.26) 


V_( k , a ) IAI-^aa) 


lAI-jCaa) 


X_(k) 2ir(a - k)x_(k) 2ir(a - k)x_(a) 

* (u7IT) 2 ((u7IYf - - 1K ^ (k ’ 


a - 


U(a) 


+(uu)) 


a) ~1AI,(aa) 

— + 


[co - kU( a) ] 


J + (k) 


2 // \ 2 
u 1 



2 tt ( ot - k)x_(a) 

\ 

2 


_ vk 2 fk 2 - 1/a 2 ) 1A ¥* a) T U»aW UiiTaW ' 1/a 

(a - k) [u - kU(a)k + (k)2ir («- 5c + (ufaj) [ " - ktJ < a )] 2w 

- = E(k) (4.27) 

2irJ + (k) 


where E(k) Is an entire function which, by the same arguments used In 
Chapter 3, Is a polynomial. 

Replacing E(k) with a polynomial and returning to the split 
functions x + by equation (4.21) we can write the solution to the 
Wlener-Hopf equation as 
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V_ ( k , a ) = (k - a.) x_( k ) 


[• 


b k 


+ . . . + q 


■4 


(aa) 


- k) 


1Al l( aa)x'(a)x_(k) | 1AI l( aa - 1/a 2 ) (k - a )x_(k) 
2irx + (a) + / w x2 


x + (u“a))(“ • u“a>) 2ir 


[to - kll(a) ] 

(4.28) 


and 


P + (k,a) = 


Kk - 


a ) 5 C + ( k ) L. a "| AI (aa)x (ot)x (k) 

; a + b k + . . . ♦ q k q - — 1 i — 

L J ~ \ 


2irx, (<*) k 

T 


AIi( a a )Y k 2 (k 2 - Va 2 ) AI l( aa) Y (^) 2 ((^) 2 - !/a 2 ) (k - a)Xt (k) 


2irk^(a - k)[a> - kU(a) ] 


* + fey)(“ - ufil) 2ir t“ - kU( a ) ]k‘ 


Y C (0) (4.29) 

2ir 

As In Chapter 3 we will choose the coefficients a', b', etc. so 

a 

that (1) P + Is analytic In the upper half plane, (2) there are no 
poles at k = 0, and (3) we obtain the least singular solution near 
x = 0. 

A 

The singular part of P near k = Q can be found to be 


1 


AI',(aa)x (a)x.(0) ACax.(O) 
- 1ax + (0)a - — 1 + 


la (- ax + (0) + x + (0)) - 


2irx + (a) 

Al^aaJx (a) X+ (0) 


1<o 




- 1ax + (0)b + 


i i 


2*X + (<*) 


AC / ' \ 

— L(-ax + (0) + X + (0))- ax + (0)U(a) 

to 1 x / 


(4.30) 
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In order to eliminate the pole at k = 0 we must set 


a = 


‘A 1 1 (®a)x (<*) ACa 


2irx + (oi) 


I 0) 


(4.31) 


and 


s' .=i 


AC ACaU(a) 
la - t~ + 


IG) 


I 0) 


(4.32) 


where 


C = 


U l (aa)Y WaW VVU(a) 


- l/a 1 


(4.33) 


2-irx, 


(O 


+ \U(a) 


Ot - 


U(a) 


This choice for a and b ensures that the pressure decays 

* 

at upstream Infinity and that p + Is analytic In the upper half plane. 
With these conditions satisfied then we can obtain the least singular 
solution by setting the remaining coefficients equal to zero. Hence, we 
can rewrite equations (4.28) and (4.29) as 

, , 1 A I , ( ad ) 1AI,(aa)x (oOx_(k) 

VJk.a) = (k - a) x _(k)(a + bk) 


2ir(a - k) 


2irx + (a) 


AC(k - a)x_(k) 
[co - kU(a) ] 


(4.34) 


and 


P + (k,a) = 


1(k - a) x+ (k) 


2 . . 2 


, X Al^otaJx (a)x + (k) 
a + b k ] - ~ 

2irx + (tt)k 


Al^aajyk (k - l/a") +1AC(k - d) x+ (k) yC ' {Q) 

2 , . V r . . V + r v n. 2 2ir 


(4.35) 


2 irk (a - k)[co - kU(a) ] [co - kU(a)]k‘ 

The purpose of constructing an elgensolutlon to the problem was to 
use It to eliminate the edge singularity In the noncausal particular 
solution (eq. (3.80)). However by examining equations (3.80) and (4.34) 
we see that the least singular elgensolutlon we were able to construct 
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(consistent with the boundary conditions at upstream Infinity) has a 
higher order singularity than the particular solution equation (3.80). 

If we try to remove the singularity In the particular solution by adding 
this elgnsolutlon we will still be left with a singularity at the edge 
of the duct and hence the resulting solution will still not satisfy the 
Kutta condition (In fact It will not even exist Ip the usual sense). 
Evidently another elgnsolutlon must be added to relieve this singularity. 
4.4 The Kutta Condition 

The difference between the causal and noncausal particular solution 
must satlsy the differential equations (2.11) to (2.13) subjet to homo- 
geneous boundary conditions, that Is It must be an elgensolutlon. From 
equations (3.17), (3.78), and (3.80) we can find this elgensolutlon for 
the upwash velocity to be 


V e (x,a) 



iP 0 x_(k) 

2irx + (0) 


+ LI k^al 


e ikx dk - 


V* (a) 


a 


X + (0)I 1 (aa)x (a) 


1-j (aa)e 


1aX 


If we rearrange equation (4.34) we can write It as 


(4.36) 


v (k,a) = 





*)x (k) 


— 1 A I (aa)x'(a) 
1 


2irax (a) 
+ 


1AI 1 (aa)x , (a) x ^(k) 

2irx + (a) 


+ (k - a) X (k) 


1AI (aa) 


AC AC 

1 

+ 

— 

(0 <j! 

2ir(a - k) 


AC + Kl . 

-] AC( k - a) x _( k ) 

+ D K 

0) 

+ [« - kU(a) ] 


(4.37) 


Combining equations (4.8), (4.19), and (4.37), we can write the complete 
elgensolutlon we have constructed as 
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V (x,a) = 


i: 


-1 AI-j (aa)x' (a) 
2irx + (a) 


£ |(k - a) X _(k) - ^ + b'kj + 


x - (k )\ + V e1kXdk + Al^aa) e Ux 

AC(k - a) x (k) 


[u - kU(a) ] 


1 kx 

e 1KX dk (4.38) 


The first two terms on the right In equatlpn (4.38) are equal, to 
within a multiplicative constant, to the elgensolutlon (4.36). The 
remaining term then must by Itself be an elgensolutlon with the same 
level of singularity as the elgensolutlon we constructed In the last 
section. 

We write this elgensolutlon as 

v E ( X> a) . £ (k - „) x_(k) - [& ♦ b'k + rfu-( ' a)] ] e ,lcx dk (4.39) 


If we added the elgensolutlon (4.38) and subtracted (4.39) from the 
noncausal particular solution we could eliminate all the singularities 
by choosing the constant A correctly and we would of course arrive back 
at the causal solution. 

Naturally we cannot conclude from the discussion of this section 
that all solutions which satisfy the edge condition are causal. However, 
with the procedure used here, we were not able to construct any noncausal 
solutions which satisfied the edge condition. 

4.5 The Effect of the Mean Flow Profile on the roots of the Dispersion 
Equation 

Before closing this chapter we explore the effect of the mean flow 
profile on the roots of the dispersion equation (eq. (4.13)). The Inclu- 
sion of the parameter b allows us to trace the evolution of these roots 
as the mean profile goes from plug flow (b = 0) to a Hagen-Polsevllle 
parabolic profile (b =1). 
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For the numerical computations we write equation (4.13) In dimen- 


sionless form as 


[a - x(l - b) ] 


1 + 


xl^(x) 

I^x) 


- 2x[Q - x(l - b) ]b 


+ &I __ 0 
B 


(4.40) 


Roots corresponding to the lowest Instability modes were computed 
over a range of Strouhal and Weber numbers. Figures 4.1 to 4.5 show the 
development of these roots as b goes from 0 to 1 . In these figures 
Im H corresponds to the growth rate of a disturbance with wavenumber 
Re x\ 
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CHAPTER 5 


DISCUSSION AND SUMMARY 

We have examined the possibility of exciting capillary Instabilities 
In a circular liquid jet by forcing the flow from within the nozzle. 

This is the so-called receptivity problem. A time harmonic axial pres- 
sure gradient was Imposed on the steady, parallel flow of a jet emanating 
from a semi-infinite circular duct. 

Using a method developed In the study of plasma Instabilities we 
were able to construct a time harmonic causal solution to the forced 
problem over certain ranges of the physical parameters. In order for 
this time harmonic solution to be causal It must contain a term which 
grows exponentially In the downstream dl rectlon, In other words pn jnst£- 
blllty wave. Hence causality provides a mechanism by which externally 
Imposed disturbances can couple to Instability waves. In addition, 
causality uniquely specifiers the amplitude of the Instability wave 
relative to that of the forcing and hence the "coupling coefficient" Is 
determined. This "coupling coefficient" Is a measure of the effective- 
ness of the disturbance In generating Instability waves. The analysis 
of Chapter 3 yielded a formula for the coupling coefficient from which 
numerical values were computed for different cpmblnatlons of mean pro- 
file, Weber number, and Strouhal number. 

In the process of constructing a causal solution we found a range 
of Weber numbers (In particular Weber numbers less than around 5) for 
which a time harmonic "steady-state" solution does not exist. For these 
Weber numbers a disturbance will grow exponentially In time at every 
point In space so that the concept of spatially growing waves Is no 
longer applicable. In this case the flow Is said to be absolutely 
unstable. Since we are Interested In generating spatially growing waves 
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we have restricted our attention to the range of Weber number for which 
a “steady-state" solution exists and have not further pursued the case 
of absolute instability. 

In addition to the casual solution obtained by the method described 
in the appendix a noncausal solution to the forced problem was found. 

The noncausal solution had a delta function singularity at the trailing 
edge of the duct. In order to obtain a solution which satisfies the 
Kutta condition a constant multiple of an elgensolutlon, with the same 
level edge singularity can be added to the forced solution so that the 
singularity Is cancelled out. 

In Chapter 4 we constructed an elgensolutlon to the equations and 
boundary conditions set down In Chapter 2. However, the least singular 
elgensolutlon we were able to construct contained a higher order edge 
singularity than the noncausal particular solution. Hence the sum of the 
two still would not satisfy the Kutta condition. We were able to elimin- 
ate all the singularities by appealing back to the solutions of Chap- 
ter 3 but this Inevitably led us back to the causal solution In order 
satisfy the edge condition. That Is to say we were not able to construct 
any noncausal solution which satisfied the edge conditions. 

In constructing the elgensolutlon In Chapter 4 we derived the dis- 
persion relation for the stability problem of the doubly Infinite jet for 
the family of mean profiles given by eq. (3.21). The lowest order root 
In the fourth quadrant of the complex wave number plane was computed for 
this family of profiles over a range of Weber and Strouhal numbers 
extending the results of Keller et al. (1973) who computed this root for 
the plug flow profile (b = 0). 
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5.2 Prospects for Further Research 


The discovery that the liquid jet can support absolute Instabilities 
was an unexpected result which deserves more attention. By a detailed 
numerical Investigation of the dispersion relation It would be possible 
to pin down, for each mean profile, the range of Weber number for which 
the absolute Instabilities arise. It may even be possible to determine 
this range analytically for the plug flow profile. 

The effects of viscosity on the coupling coefficient and growth 
rates of the Instability waves have not been considered In this work and 
would make an Interesting extension of the results obtained here. 
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APPENDIX A 


A.l Introduction 

In this appendix we will outline some of the theory behind the cri- 
teria used In the main body of this worl^ to ensure that causality Is 
satisfied and Identify absolute Instabilities. The mathematics was 
worked out originally for the study of plasma Instabilities. We have 
used the results as derived by Briggs (1964) and Bers (1972) and refer- 
ence will be made to them for detailed proofs. 

The method for obtaining a causal solution Is based on the examina- 
tion of an Initial value problem where the forcing Is "turned on" at some 
Initial time. The "steady state" solution will be the long time behavior 
of the solution to this Initial value problem provided that the flow Is 
not absolutely unstable. 

A description of the method and a discussion of the Issue of abso- 
lute Instability Is given In the following two sections. Since we have 
a specific example In the present problem we will use It to Illustrate 
the Ideas. 

A. 2 Method of Solution 

We consider an Initial value problem. At time t = 0 the flow In 
the semi -Inf Ini te jet Is subjected to a pulsating axial pressure gradi- 
ent. The response of the flow to this forcing In space and time can be 
determined by an Inverse transform of the Greens function for this prob- 
lem In the frequency and wave number planes as 

, t . f f rr n 1kx 1wt hl- 

g(x,t) = J J G(w,k) e e dk do> (Al) 

L F 

or as a function of frequency as 

g( x ,co) = j G(w,k) e 1kx dk (A2) 

F 
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This Greens' function win have poles at the zeros of the dispersion 
relation In the k and w planes. The Inversion contour In the fre- 
quency plane must be a line such that no zeros of the dispersion relation 
(complex <o for real k) exist above It, since only then will the con- 
dition of causality be satisfied (for t < 0 the Integral Is evaluated 
by closing the contour In the upper half <o plane). The Integration con- 
tours In the k and w planes are shown In figure A.l. Therefore, we 
need to obtain the solution in the frequency plane along this contour L, 
that Is for the frequency having a sufficiently large positive Imaginary 
part. 

We are interested In the asymptotic response of the flow a long time 
after the forcing was Initiated. The asymptotic response can be obtained 
by moving the contour In the co plane as close as possible to the real 
axis (since, In the absence of absolute Instabilities the behavior of the 
solution for large t Is determined by the pole of g nearest the real 
axis In the <o plane). If we can move the contour all the way to the 
real axis, that Is If we can analytically continue g to the real axis, 
then the long time response will be due to ^ pqle on the real axis which 
corresponds to the frequency of the forcing. In this case then a time 
harmonic, "steady state" solution exists. 

Now the poles In the k plane are related to to through the dis- 
persion relation, say A(k,co) = 0. As we let to approach the real 
to axis it may happen that certain poles In the k plane cross the 
real k axis (In fact. If Instability waves are to be generated by the 
forcing this must happen). If this happened g(x,to) would not be ana- 
lytic unless the contour In the k plane Is deformed so as to continue 
to Include (or exclude) any poles which might be Inclined to cross the 
real axis as Im co -> 0. This Is Illustrated In figure A. 2. 
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The spatial response of the flow Is obtained by closing the Integra- 
tion contour In the upper (lower) half k plane for x > (<) 0. Any 
poles which crossed Into the lower (upper) half plane as Im w -> 0 (and 
which remain Inside the contour due to the deformation) correspond to 
growing Instability waves In the "steady-state" response (see eq. (A. 2)). 

In short then, to obtain a causal solution, we need to obtain the 
transform of the solution to the Initial value problem with the frequency 
having a large positive Imaginary part. This solution will appear form- 
ally the same as If we were seeking a time harmonic solution from the 
start. Hence we can proceed from the time harmonic form while keeping 
In mind that the frequency has a large positive Imaginary part. We then 
let Im to -> 0 and keep track of any poles In the k plane which cross 
the real axis and In this way obtain a causal solution for to real. 

A. 3 Absolute Instabilities 

It may happen that the analytic continuation just described cannot 
be carried out all the way to the real w ax1 x s. One way this could 
occur Is If two roots In the k plane, one originating In the upper and 
one In the lower half plane, merge for some complex with Im w > 0. 

As can be seen from figure A. 3 the Integration contour In the k plane 
will become "pinched" between the two poles and the deformation on which 
the analytic continuation rests cannot be performed. In these situations 
(as shown by Briggs and Bers) the asymptotic time response Is dominated 
by a term which grows exponentially In time at every position In space 
and the flow Is absolutely unstable. 

When these absolute Instabilities are present there Is no "steady- 
state" and the concepts of spatial Instability break down. 
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Figure 2. 1- Flow geometry. 
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Figure 3. 1- Movement of the roots of the Kernel function as 
Imfivaries. 
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Figure 3. 3- Locus of complex root of x for b = 0, (r = 10. 



Figure 3. 5- Locus of complex root of x for b = 0, (3 2 = 30. 
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Figure 3. 7 - Locus of complex root of % for b = 0, = 50. 
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Figure 3. 9 - Locus of complex root of x for b = 0. 2, fr = 20. 




Figure 3. 11 - Locus of complex root of x for b = 0. 2, p 2 = 30. 




Figure 3. 13 - Locus of complex root of x for b = 0. 2, f} 2 = 40. 







Figure 3. 20 - Movement of roots of x In the case of absolute 
instability. 







Figure 4 1 - Lowest order root of the dispersion relation, b = 0. 



Figure 4 2 - Lowest order root of the dispersion relation, b = 0. 2. 
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